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ABSTRACT 


A finite element method was applied to the unsteady 
transonic small disturbance equation and integrated until 
the solution converged to the steady state for a thin non- 
lifting airfoil. The method of weighted residuals was used 
to formulate the finite element equations, and Houbolt's 
method of central differencing in time was used to integrate 
these equations. 

A secondary investigation applied the steady transonic 


small disturbance equations to a converging-diverging nozzle. 
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I. INTRODUCTION 


Transonic inviscid flows past a smooth airfoil may be 
expressed in terms of the velocity potential ® satisfying 


the transonic small disturbance equation, 


CS ee ee (1) 


fursecquatiren presents two mayor ’difficultiesf.1) it-is non- 
linear and 2) it is of the mixed hyperbolic-elliptic type. 
Analytical solutions to non-linear equations are difficult 
to obtain. One must normally resort to numerical methods. 


2 2 


When the coefficient (l - M - Mey 1)9,) in equation l 


is negative, the flow is supersonic and the equation is 
called hyperbolic; otherwise the flow is subsonic and the 
equation vis Giliptic?: > vikesforms Sef *theetwoysolutions ‘are 
fundamentally different. Hyperbolic equations admit both 
discontinuities, which propagate only in characteristic 
directions, and the presence of shock fronts. Elliptic 
equations, on the other hand, require that the dependent 
variables and their derivatives be continuous and that a 
Ehange ti any part of the flow field affects every other 
part. Many non-linear elliptie equations are solved by 
apprepriate relaxation iteration techniques by. casting the 
equation in Poisson's form with the non-linearity acting as 


EMe a@tivine force. Solutions to hyperbolic equations are 


usually obtained by the method of characteristics or by 
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finite difference marching techniques. which use an. artificial 
WiScOSity to sepresent, the average jump conditions across the 
shock wave. In mixed supersonic and subsonic flows, normal 
numerical procedures break down because the boundary between 
the two regions is not known a priori. 

Finite element numerical techniques have evolved as a 
powerful tool in obtaining approximate solutions to a wide 
variety of engineering problems, particularly ones with 
Neumann-Dirichlet boundary conditions, i.e., elliptical prob- 
lems. They offer several outstanding advantages. Some of 
chese=are: 


1) Non-homogeneous problems may be treated with relative 
Simplierty. 


2) Complex geometries may be modeled with relative ease 
Since the elements can be graded in size and shape to 
follow boundaries of arbitrary shape. 

3) Once the finite element model has been established, 

a variety of problems can be solved by supplying the 
computer with the appropriate data. 

Chan and Brashears [Ref. 5] developed a finite element 
computer program for steady transonic flow over a non-lifting 
airfoil. This program uses the least squares method of 
weighted residuals to approximate equation 1 by a system of 
algebraic equations, and assembles the equations in a special 
May te decount for the hyperbolic region of flow. This tech- 
nique prevents the influence of downwind nodes from propagating 
upstream in the supersonic region. 

ie euepese Of Chys thesis 1s to investigate the possi- 


bility. of speeding the convergence of a solution by transform- 


ing the steady transonic equation to the unsteady equation 
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ana imtegrating over time until the time dependent terms 
vanish. and to extend.the program of Ref. [5] to the tran- 


sonic region of a converging-diverging nozzle. 
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1 DISCUSSION OF THE FINITE ELEMENT APPROACH 


In a continuum problem of any dimension, the field 
Variable, whether it is velocity potential, velocity, temper- 
ature, displacement or some other quantity, possesses infin- 
itely many values because it is a function of each generic 
Home an ehe Solution region... Consequently, the problem. is 
one of an infinite number of unknowns. The finite element 
approach subdivides the solution domain into a finite number 
of subdomains called elements and expresses the unknown field 
variable in terms of assumed approximating functions within 
each element. The approximating functions are sometimes 
eaided anterpolating functions and are-defined. in, terms of 
the values of the field variables at specified points called 
nmeadges Of nodal points... Nodes usually, lie on, jthe element 
boundaries where adjacent elements are considered to be con- 
Hected. In addition to boundary nodes. an,.eclement may also 
contain interior nodes. The nodal values of the field variable 
Sonehe Sc erNOlating, function forthe elements completely 
define the bahavior of the field variable within the elements. 
Once these unknowns are found, the interpolating functions 
define the field variable throughout the assemblage of the 
element. 

Clearly, the nature of the solution and the degree of 
accuracy of the approximation depends not only on the number 


and size of the elements used but also on the interpolating 
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functions selected. Interpolating functions may not be 
chosen arbitrarily because certain compatibility conditions 
must be, satisfied... Qften such functions are selected: so 
that the field variable or its derivatives are continuous 
across adjoining element boundaries. Once the problem is. 
formulated in terms of individual elements, the contributions 
of each element may be assembled to define the entire solu- 
tion domain. This means, aon example that af we: arestrea ting 
a problem in stress analysis we can find the force-displace- 
ment or stiffness characteristics of each element and then 
assemble the elements to determine the stiffness of the whole 
structure. Finite element solutions are not, of course, re- 
Senterea to structures problems, but the matrix of equations 
defined by the interpolating functions and the nodal field 
Warrapves is stills reterred to as the stiffness matrix 
regardless of the field variable in the problem. 

SolLgcrons co CcOMeinueus problems by the finite element 
approach always follow a systematic step-by-step process. 
ines process is ‘completely general to’ the finite element 


method and it is outlined below. [Ref. 4] 


L, Drseretize the continuum. 


The first step 2s to dimide the: solution domain into 
elements. A variety of element shapes may be used and one 
or more different element shapes may be used in the same 
region. The type and number of the elements used in a given 


problem are a matter of engineering judgement. 
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we select the-interpolating functions. 


The next step is to choose the type of interpolating 
function to represent the variation of the field variable 
over the element. The field variable may be a scalar, a 
wector, or a higher order tensor. Often polynomials are 
selected as interpolating functions for they are easy to 
integrate and differentiate. The degree of polynomial 
chosen depends on the number of nodes and the nature and 
number of the unknowns and the continuity requirements 
imposed at the nodes and the element boundaries. The un- 
known quantities at the nodes may be assigned to the field 


WVaeriaple and their derivatives. 


3. Find the element properties. 


After the elements and their interpolating functions 
have been selected, the matrix of algebraic equations which 
express the properties of the individual elements must be 
determined. Several methods are available for this task. 
iiese are: the direct approach, the variational approach, 
the method of weighted residuals and the energy balance 
approach. Reference [4] is a good source of information 


on the various techniques. 


4. Assemble the element properties to obtain the system 


equations. 


ieetina the properties of the over-all. system, the matrix 


equations expressing the behavior of each element must be 
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added to the matrix equation of all other elements. The 
basis for this assembly procedure stems from the fact that 
connecting elements have common nodes and the field variable 
must be the same for each element sharing that node. 

Ae this “point the boundary ‘condat ions: *forsthe system of 
equations must be accounted for and the system of equations 


must be modified before it is ready for solution. 


a., solve the system of equations. 


The assembly process of step 4 produces a set of simul- 
taneous equations which can be solved to obtain the unknown 
field variables. Linear equations have a number of standard 
solution techniques readily available, but solutions to non- 


linear equations are more difficult to obtain. 


6. Make additional computations if desired. 


Important parameters, such as pressure coefficient in 
aerodynamic problems, may now be calculated from the values 


of the field variables. 
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TIt. THEORY AND BASIC EQUATIONS 


A. 2 STEADY «TRANSONIC. FLOW 

Chan et al. [Ref. 5] developed an algorithm to analyze 
steady transonic flow over non-lifting thin airfoils. Bound- 
ary layer effects were ignored and the imbedded shock wave 
was assumed to be weak. These assumptions are consistent 
With transonic small disturbance theory which can be expressed 
mathematically by the following expressions. 


(l - M2 - MECy + 1)0,)0,y * byy 7 0 a) 


Boundary conditions - 


V- o-= 0 ae Antana ty (2) 

v = (1 + uj)dg/dx on. the airfioit (253) 

u = 0 on the line of symmetry 
(4) 


where g(x,y) defines the airfoil and dg/dx describes the 
airfoil slope. 

_ The above expressions appear in their dimensionless form 
where > = perturbed velocity potential and the perturbed 


velocity components in the x and yY directions are respectively 


defined as 


u = a 
v= % 
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M. = freestream Mach number and ~ = ratio of specific heats 
Witter fOredirers taken tetbe By4. The physical coordinates 
x* and y' and the velocity potential %' are related to the 
dimensionless quantities by 
Sec, yee the; a Ole 

mMnere © 15 Ehe Chordvength of the airfoil and Uy is the free- 
Stream velocity. 

Once the flowfield solution is determined in terms of the 


perturbed velocities, the secondary unknowns are subsequently 


calculated. These include: 


: Oe ee 
a= (Feu -v) + G12 (5) 
ae 
M = = (6) 
Epes shee Lonel EMned Ot 1) 
‘9 [2 a Ghe] y/(y +1) 
pe eae? 
eet ee) ee | a (8) 
co U U co 


Pm chic. above, Uo — 1 is the normalized freestream velocity, 
a= local sound speed, p = local static pressure, M = local 
Mach number, V = total velocity ; Py = stagnation pressure 


and C5 = Pressure cockriezent. 


B. UNSTEADY TRANSONIC FLOW 
Unsteady transonic inviscid flow may be expressed in terms 
ge tne wMelocity potential O(x,y,t).to a first order approxima- 
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which has the same non-linear coefficient retained for steady 
transonic flow in Equation l. 
Boundary conditions require that the disturbances vanish 


far from the airfoil, 


6 = 0 eo. = 0 at infinity (10) 


and that the flow remain attached to the body. Let B(x,y,t) =0 
define the body at any instant. The surface tangency restraint 
may now be expressed by the substantial derivitive DB/DT 


vanishing. 


DB/DT = B. #OClest > SBE & o By (Gals) 


If the body remains stationary Be = 0, and the tangency 
condition becomes the same as in steady flow, namely 

MY = (1 > w)ide/dx (Zz) 
where dg/dx represents the airfoil slope. 


The pressure coefficient for isentropic unsteady com- 
pressible flow is defined by 
= 2 Dole vee 2 7a) AL) 8 
epamdasor {01 =o Me 2h 2G th + OD] 1 
ee (13) 
Expanding by the binomial expansion and retaining only the 


first order terms gives, 


Co c0,. >) 24 (14) 


Pp ie 


ay 


_ 7 4" J 


; . 

2 ayy Son _ 4 ee ae 

= ++! M a Me “ee aH hye © f rns — 
‘ <i) 5 


jay Ineisitieus ~cagne t “nok ‘sie OE ee Bei 
: 

| soPeus@ af eel? ofae ene 

nS, 

sid 3587 Sehupée candies Sonia 

i = 


tzotste od3 mor 1a) 
ia is 


ts 


} fi 


23 bedoette diames wold eft goede bas 


i ie s\T .3esteay vee ge vbood ods. enit b 
es 
cijoasecue of 46 beezetqxe sd won Vee 
» pabdz a 
~ ome ea 
* = t. % 2 ° 7” 7 Ta \ aq 


; : - 7 
7 Soe .0 = | ytameisete enbemes ybod ede Zi! 
pal - a 
51% koje ci *8 SORE Sad. gau6sed aos ooD 
7 : vo — 

eh\gbiaw «+ f) = y 


: 7 
-sqoule sfektss ene aseebutyey xb\gb roi 
7 


mu 2iqorImez) tot saelosetee etwee org on 


“a benlisb * 22 wolt otai 2297 
; ~) : ‘gaits + _ 
‘9+ gtey s) yu a a) ee wate? 
v 


is a 
co ne Moa pets 
7 A veg am 2% 7 sto te" 


i ie a 
® L” 


7 


a 
: - i 


IV. FINITE ELEMENT FORMULATION 


The method of weighted residuals is a technique for 
approximating solutions to linear or non-linear partial dif- 
ferential equations and it is the basis for the finite element 
formulation of the transonic small disturbance equation 
(Equation 1). This procedure involves assuming the general 
functional behavior of the field variable which would approx- 
imately satisfy the basic equation and boundary conditions. 
suDStituting this approximation into the original differential 
equation results in some error called a residual. A system 
of algebraic equations results when a weighted average of the 
residual 1s forced (to vanish as it as averaged over the entire 


domain. 


A. STEADY FLOW 


The approximate solution to equation 1 is assumed to be 


> = N.d; (15) 


alert 
where .N; are the interpolating functions which exhibits the 
behavior of equation 1 and >; are the undetermined parameters 
at the nodal points. 

When oe SuUpStituced anco equation LT, the resulting 
residual is 


ee [tM - ME CY FLING IN sex NIA (16) 
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The weighted average is determined by multiplying the 
residual R by m linearly independent weighting functions W; 
and integrating over the elemental domain. Forcing this 


residual to vanish yields, 
ffW;RdA = 0 = 1b £6 m 


Chan et al. [Ref. 5] found that when the weighting function 
Ws Pon pequation. 1 15 chosen to be 9R/ 36; the Tesulting matrix 
is positive definite and well conditioned. This choice of 
weighting functions is referred to as the method of least 
Squares because it is equivalent to minimizing the square of 
the residuals summed over the domain with respect to the un- 


determined parameters. That is, 
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axy/a % = ssaR/ao;RdA = 0 (17) 


Integrating over the domain produces the system of 
algebraic equations 


= 0 (18) 
where the elemental matrix oe is defined as 


95, = FIQTES dA (19) 


With Q; and P. equal to 
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B. UNSTEADY FLOW 

Development of the unsteady flow finite element equations 
is similar to the procedure used to formulate the finite ele- 
ment equations for steady transonic flow. The least squares 
method of weighted residuals is again used but » is now a 
function of time as well as the spatial coordinates x and y. 


The approximate solution has the form, 


d = Ni (x,y)o, (t) (20) 


substituting $ in the unsteady transonic small disturbance 


equation, the weighted residual becomes, 


a= SiR, + R, + Rz) “dA (2) 
where 
Be tM Me Gedo N IN Nee eas 
at 00 co SY Ok | ere. Sear Soe Ie5ye 
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Zz se 
Re = -M'N.»&. 
3 jj 


Expanding equation 21 yields 
“ 2 2 2 
ieee Ly Be aR. + 2R,R, + 2R,Rz]dA €2Z)) 
and minimizing with respect to the undetermined parameters 
5 the following system of algebraic equations result, 


b; = 0 = 2ss9R;/99; (Ry + Ry + Rg]GA 


aRy/ 995 ae ge 


20 


~ 
P| 


vearmTs 7s bay saz ‘os —_ 


_— 
teluwrol az boars Srebeong — miimt 
ol? Sincanery “yonoae 9 sups : 


iS cs ng $95 @L slaubbews | yes aes: bods 
ite 


inftisqe ed? 4s itew a a0 33 — seks: >t 
tot sds eed roisiter etnatxe — 
7 _— 

(9) 54 cea) i= i 


- 
A - 
Mie ILROZtsTs yaar eee on3 at } as jised & 
ied 
,eemoned Snub _ tendater Sued mot 


one 
se Ss 
ces 


be” +A ae = y ay 
: y, 


\ 


ee 
oe) 


_ ig GPn 2 a i+ a 


ehlery is, nassaine 


a 
Ro + AAS Saw tee ie ay 


. 


ru 


nol tees isndeade 7 


Ao | ls oF “ae 


where P. has been previously defined in the steady finite 


element formulation 


The above equation may be rewritten in the form 
SweGe 2 Se ole SMe 6S 0 23 
ij°3 ij°3 ij°3 cP 


The stiffness matrix Si; is the same as that developed 
for the steady transonic equation and the damping Ce and 


mass Se matrices are defined below. 
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V. ELEMENT DESCRIPTION AND ASSEMBLY OF EQUATIONS 


A. ELEMENT DESCRIPTION 

The basic element used in the finite element program is 
the non-conforming cubic triangular element developed by 
Bazeley etal. [Ret.. Z|). Also used in the program are the 
quadrilateral and trapezoidal elements constructed from these 
triangular elements. These three types of elements can be 
mixed and used freely in the entire flow region except that 
only trapezoids should be used to cover the supersonic and 
mixed region in order to enact the special assembly procedures 
required by the hyperbolic equation which describes the flow 
im that region. 

The basic triangular element is shown in Fig. 1, which at 
each vertex has the velocity potential and the velocity com- 
ponents as the undetermined parameters. This type of element 
was chosen because both Dirichlet and Neumann boundary condi- 
tions can be imposed with equal convenience. In addition, 
because velocity components are used as primary unknowns 
secondary parameters, such as Mach number and pressure coef- 
ficient can be calculated directly without resorting to 
numerical differentiation, which would produce additional 
errors. 


In the element, the approximate solution is assumed as 
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Pagure 1 — Triangular Element 


ZS 


In which 9,;'s are the nine undetermined parameters of » and 
N; are the znterpolation functions which are expressed in 
terms of the area coordinates. 

The shape or interpolating functions and their first and 
second derivitives are defined below. 
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Quadrilateral and trapezoidal elements as shown in Fig. 2 
are also used in the present program, the former is used in 
the subsonic region and the latter in the mixed and super- 
sonic fegion. The element matrix for the quadrilateral ele- 
ment is obtained by combining appropriately the matrices for 
two triangles, while the matrix for trapezoidal element is 
obtained by averaging contributions from two left-running and 
two right-running triangles. The averaging process is designed 


to remove the bias. effects inherent in the quadrilaterals used. 


B. ASSEMBLY OF EQUATIONS 

Straightforward application of the finite element assembly 
technique to transonic flows would fail (the solution diverges) 
because this would allow disturbances to propagate upwind in 
the supersonic region of flow where the governing equation is 
hyperbolic. Hyperbolic equations have a time-like dependency 
in that the solution at the downwind station is affected by 
the upwind station but not vice-versa. Assembly techniques 
£0r 2 transonic flow finite element program must take into 
account this time-like dependency. If the x-axis is taken 
aoa cime-like direction in the supersonic region, the ele- 
ment matrix may be assembled in a way similar to a backward 
finite difference operator, which has been successful in 
solving hyperbolic equations. 

Consider the rectangular element sketched below with the 
upwind station I and the downwind station II, each having two 
nodal points with the element type chosen. The element 


matrices can be constructed in the usual manner. . 
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However, before assembling the element matrix into the system 


matrix the non-linear coefficient of equation 1 is evaluated. 
Gaol Me i ty + 2) u 


the sagen of the. coefficient being positive, zero, or 
Negative detines the equation as elliptic, parabolic, or 
byperpelic. if € is non-positive for all nodes in the element, 
the rows representing the improper-downwind influence on the 
solution at an upwind station are ignored during assembly. 
This feature is automatically applied in the program requiring 
aimee !ttele Care in-arrangement of the nodes of the element. 
In the anticipated supersonic region, element node points 
should be arranged in the order as indicated in figure 3, 
Starting with the upper left. corner and proceeding in the 
eoumter-clockwise direction. in the elliptic région, i.e., 
where the coefficient is positive, no special assembly tech- 


nique is invoked. 


Ga ETERATIVE PROCEDURES 
With the equations assembled and the proper boundary con- 
ditions imposed, the system of non-linear algebraic equations 
is solved by iterative procedures in the form 
ys UR) © 2 
S45 (9) 9; Ls (23) 
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£0 solve for the solution in the n Lperation. The fLunction 


o is defined as 
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Figure 3 - Nodal Arrangement for Supersonic Region 
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in which the under-relaxation factor 9 is in the range 
ea. er subsonic flow 9 = 1, which as simply a succes - 
sive approximation, yields good results, but it is necessary 
to under-relax somewhat with 9 approximately .5 for super- 
errcreal flow. Generally, a smaller relaxation factor will 
Make ehe solution more stable but it will tend to slow down 
the rate of convergence. 

Equation 23 is subject to the convergence criterion that 
ene, change in local Mach number between two, successive iter- 
aetons 21s less than a prescribed’ value) «& at alle nodes: an: the 
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VI. INTEGRATION OF UNSTEADY FINITE ELEMENT EQUATIONS 


The unsteady transonic small disturbance equation (equation 
9) when suitably reduced to a finite element approximation 


appears in the forn, 


Gros ve SC. G+ SM. .6. 2 B | 26 
ij? j ij° 3 ij? J E26) 


This equation is analogous to a damped spring mass system, 


hence? ~S:= , a: : are respectively referred to as 


the stiffness, damping, and mass matrices. 


and SM. 
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Mathematically, equation 26 represents a system of second 
erder differential equations with constant coefficients, which 
ean we solved by standard numerical procedures for differential 
equations, such as Runge-Kutta or Milne methods. However, 
this would be a very costly technique if the coefficient 
Matraces dre. very) large. in practical finite element, analysis 
there are a few effective methods which take advantage of the 
banded matrices usually encountered in finite elements. One 
such method is the direct numerical integration method. 

Direct integration involves a numerical step-by-step 
procedure aimed at satisfying equation 26 only at discrete 
tmie intervals At apart.and not over all-time t. Conceptually, 
direct integration is a finite element method in space and a 
finite difference method in time. Examples of direct inte- 
gration are the central difference method, Houbolt integration 


and the Wilson method. The first two schemes are finite 
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difference schemes whereas the latter is a linear accelera- 
tion method. Linear acceleration integration assumes a 
Hinear yariation of acceleration from time t to time t + At. 

Central differencing can be very effective in the solu- 
tion of many dynamic problems especially those that involve 
a large system of equations. However, this method is unstable 
Por ali time steps larger than a critical time step: 

Houbolt integration is an implicit finite differencing 
method related to central differencing, only it has the ad- 
Vantage of being stable for all TIME STEPS. 

The Houbolt method was used to integrate the unsteady 
finite element equations because of this stability. 

Houbolt integration uses the following finite difference 


expansions: 
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which are two backward-difference formulas with errors of 
order (at) ?. 

The solution of Oi teat 
time t+At equation 26 becomes 


Must Satisfy equation 26 and at 
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Substituting the finite difference formulas for 6 ae: and 
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rearranging all known vectors on the right hand side, the 


solution for us is obtained, namely: 
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Wkere the’ constant integration coefficients ‘are: 
i 2/at? 
a, = Hy/ OAC 
ay = Sy iAge 
a, == S/he 
ay = 2a 
ac = -a,/2 
a, ay /2 
as = az/9 
Equation 27 may be written as 
SEG ae = RE; (28) 


where the effective stiffness matrix ra and the effective 


load vector see are defined as: 
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Shee = (Sea i “2 


1j LJ see i ca ee 


ee ee et ceene © Maj tocar? 


+ 2 a p ? 


ee GS eee 87 ft ne 


Accurate knowledge of the vectors $. and ¢. 


ee ‘j,t-2At 


meetedudtred tO yield -an accurate, solution for 2, t+at and 
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normally the Houbolt integration scheme requires a special 


Seameiue procedure to determine the initial two vectors 


eae tt © s 2K 


this problem is to integrate the equations until they con- 


However, since the primary interest of 


verge to a steady state solution, it is not necessary to 
obtain \an.accurate time history of the flow. Errors induced 
by the inaccurate starting vectors will vanish as time 
approaches infinity. ,[Therefore, the starting vectors may be 


chosen somewhat arbitrarily. 
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VII. CONVERGING-DIVERGING NOZZLE 


S=tF:y»Shen- [Ref.-10] demonstrated the feasibility, of 
calculating compressible flows through a converging-diverging 
(Laval) nozzle by dividing the region of calculations into 
Ehree patches, a subsonic region, a supersonic region and a 
transonic one, of course bounded by the other two regions. 
The locations of the boundaries for each region were chosen 
arbitrarily provided the sonic line is bracketed by the 
subsonic and supersonic boundaries. 

Two different finite element formulations were used for 
the subsonic and supersonic regions, but Shen [10] resorted 
PMditaiveical approximations to cover the transonic patch. 
Rhitscenestricced, the calculations to nozzles with small throat 
curvatures because no analytic solutions exist for nozzles 
wrtn laree throat curvatures. [ft is conceivable that STRANL-II 
could be adapted to provide a continuous solution throughout 
all three regions. 

Outside the transonic region of flow the governing small 


disturbance equation is 


2 is 
at (29) 


This holds for both subsonic and supersonic flow. Comparing 
equation 29 with equation 1, the transonic small disturbance 
equation, we notice that only the non-linear coefficient 


M2 (y + 1)u distinguishes the two equations from each other. 
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This coefficient becomes negligible when the Mach number 
becomes “ess jthan..8 or greater. than 1.2..- With this. con- 
sideration in mind, it was assumed that equation 1 would 
adequately describe the flow through the nozzle and that 
the finite element formulation developed for the non-lifting 


airfoil would apply to the Laval nozzle. 


A. BOUNDARY CONDITIONS 

Two solutions are possible for a converging-diverging 
nozzle: 1) Symmetric flow, where the flow is subsonic through 
the domain, except for a small supersonic region near the wall 
in the throat, and 2) Asymmetric solution, where the flow 
accelerates to sonic velocity in the. throat and. then continues 
to accelerate to supersonic velocity in the diverging section. 
Different boundary conditions apply for the two solutions. 
Port the symmetric case, both tnlet and exit velocities must 
Bemsapcetraced. Inlet and exit velocities are Equivalent in 
fie Supsonic solution.~ The supersonic solution requires that 
Omiy the anlet, velocities. be specafied. If the exit veloci- 
ties are also applied, the problem is overspecified and the 

/ 

solution may not ‘converge. 

Velocities at the inlet and exit are not uniform in the 
ye @areetion, therefore the disturbances cannot be set to zero 
asetm the case, of, the non-liftang/airfoil.. Boundary velocities 
must be calculated by solving equation 29 analytically. 

Equation 29 is a linear equation which can be mapped to 


Laplace's equation, 
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By tetting y' =./ (1 -M_) y. Laplace"s equation is easily 


solved for the case of the hyperbolic nozzle by transforming 
from cartesian coordinates to sibaie coordinates. This 
transformation simplifies the solution because the stream 
lines must be hyperbolas to follow the nozzle boundary and 
therefore follow the hyperbolic coordinate y = constant. 

Pe ee vellaperc, coordmates 4 “and y are’ ‘chosen such. that 
tne curves w = constant are ellipses andthe v curves are 
hyperbolas, then the velocity potential which satisfies 


bapiaee's equation for-a hyperbolic nozzle is simply 
od = Au 

where A..is a constant of integration. The stream function is 
y = Av 


The transformation w = u + iv = cosh +(2z/a) gives rise 


to the elliptic coordinates 


ya i 2 awecosh cos MU, x. = b/ 2 a Sinh wesin Vv 
2, 2 

Ge L/2 a [cosh t— sin: | 

a asf Cy * ADS + x? 

ge (y - ae + x? 


Solving for yu and v produces 
~ sul 
u = cosh Pury + r,)/a] 


y= cos’! Ken) al 


The nozzle boundary is defined by Me = constant, which 


along with the equation for the nozzle wall in cartesian 


5 


vitess 2 aobasups 2% ad 
poimoteres3 wf efzsom 2 
etd? ,@8*srfheneo an at ft 
wnerte sult ventebell iii gs 
fine yteinued’ aLeron adit welled = son - 
‘tnaténos = y saimERVeds pebeegalecte: Sal 
tad? afaue oe rodo eth. y ined  2SamREbT OND: ak 
ete zeviua ¥ oH?) ine poayetia: ets julieys 
seiteimee dota Ls iatissog vrtoulsy $82) aero 
yigqnize 24 sitzes. otiodgeqes © +? ne au fe 


we =" sé . id 


™ 


i noétonu? sieette- sat .noltergsanh io sasteads & aioe 
; i ~ - - 7 


» U : : "yy 
Poe 0 : : pi 


7 
y 


a iy a 
: - _ fe = Y - 7 
eeis eeavig Le\2S)” dees = Va > aie HOLT amiodel e132 of 
< - ‘ ; a5 : ; 
a _—. “ 


vA. > y 


i _ ; 4, 
Ujiiw , ca t2n0S a bl % 


malexrias,.cl Liaw § 


coordinates, y" = £(x), implicitly defines the constant a. 
Substituting for yp in the velocity potential produces, 


1 


= cosh [(r, + r,)/a] 


from which the velocities may be determined. 
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The constant of integration A may be determined by speci- 
fying the flow rate through the nozzle, but when the inlet 
Welocities are normalized with respect to the freestream 
velocity (U,) A is factored out of the problem. 

Velocities for compressible flow can be solved by mapping 
paee te the physical coordinate. system (x,y plane). 

Other boundary conditions are universal to both problems. 
These are: 

@= 0 on the line of symmetry 


(1 + w)dt£/dx at the nozzle wall 
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F(x) defines the nozzle boundary in terms of 2 ratio of the 
throat semi-height as a function of x. The throat semi-height 
is taken to be 1 for convenience. 

Pressure ratio, sound speed, and Mach number are calculated 


as before by, equations 5 through 8. 


59 


. : , os 
a : i 7 : 
7 = , ., oa = ? - ry ae 
wala ite civo2r = Fe amiey: is cmhnind oft 
: ; 7 
+ 36 gated 
ee eo Se. 
6 . evno> Pit , oe 


: coy i 
- . ae : 
oils a. “= oa = VF ie 
beselesia; ese tredmen dog bre Be ne bgt yoke 7 
: “a 7 7 “ - 


oe des E atio t tm ipo vG ® 
ar : Fi 7 ve 7 a 


meni ¢ 


_s re 
 tdgted-imee 2eo7ds ad? 


VIER. DISCUSSION OF RESULTS 


A. TIME INTEGRATION TO A STEADY STATE SOLUTION 

As stated before, the Houbolt method of integration is 
Sudke for adi timetsteps. "Results of the °test*cases bear 
this out with the larger time steps providing the most rapid 
convergence to a steady state solution. Time steps were 
Eepieecerom t = .1 to t = 100. Time steps larger than t = 100 
were not attempted because as t becomes too large the influ- 
ence that the damping and mass matrices have on the effective 
stiffness matrix becomes negligible compared to the stiffness 


macri., That cis: 
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The starting solutions were chosen somewhat arbitrarily. 


o was chosen to ‘satisfy. "the first iteration “of ‘the ‘siteady 


7o2ht 
solution 
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when the non-linear term (u) in the coefficient 


1- M2 - M@(y+1)u 
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This starting procedure proved to be superior to chosing 
Eme £irse three vectors closer to the converged solution. 


EE 6; were chosen to be the last three 


j,2at? °5,at? 0 


time steps of the previous case, the solution oscillated and 


and 6. 
eis 


eeawereed much slower than with the starting solutions: chosen 
as above. 

The stitiness, mass, and damping matrices were recalculated 
feeotecaicn time Step, using the wunder-relaxation technique 
weeseriped above. This was necessary to utilize the special 
assembly procedures invoked by STRANL-II to prevent the in- 
admissable influence of downwind nodes from propagating up- 
Seceamt an the SuUpersonie regien. 

Per Derely critical, flow (Mi= .so1) and. subsonic flow, 
an under-relaxation factor 9 = 1 (successive approximation) 
resulted in convergence to a enedng State solution after only 
three time steps. Eleven time steps were required for the 
supercritical solution to converge using the same relaxation 
Baeten-s Reducing Y fo =5 imereased the “rate of convergence 
and the solution achieved steady state after six time steps. 
Figure 4 compares the steady state solution for a 6% thick 
ceinewlar are dirtorl at Mo = -909, using “the same integration 
Methed, With the results obtained in Ref. [5]. Chan's results 


eonversed in 10 iterations ater using the results from. the 


panei, critical flow as an initial guess to the 


41 


gcieod> oF roezegie. od ar 
tofteloe As ateyano ens ox 

oe xis t2zai ot 26 60 oosnith Sia 

ona bare li ise | ieordula® ont a oe 


¢ 
neeon> 2uoitelog yaesaawe ead ditwn 


10tJaxplet~pednp). Sis unteis: « 
latoeqe ef? ocd idan 69 ae 
-ni oft tnevetg od T1-aVAgRE yd — 
-{y gtisegeroig met?  esbon ener ae 4 
»wolt cisoedee Baa (308, _ watt & 
solisemixorqgqe @Vieseooue) = & oto 
tne vette avitaiee v7at%2 \htege 4-07 
6/3 102 beteupet-o7ew’ eaete shite 
noiteaxyple: smaie saz gat eu syaevao> 0 ? 


‘ 
= 7 
o2nayrsvno 9 Ssaa ee | 


-Eqe32 omit xie sets stage 
dots: ie 302 bahar 
ijgcasiai emee ent antat 

e7iveat 2'ocdd .. [2] Aafia 


ed3 sor} 23Iueer edt 


oft 03 ve 


Supercritical solution. Figure 4 is a plot of local Mach 


numbers at boundary nodes on the airfoil. 


B. CONVERGING-DIVERGING NOZZLE 
The nozzle chosen for the test cases was the two-dimen- 


sional Oswatitsch nozzle with the boundary defined as 


ea af 2 a 


mmere the throat-at x= 2.5 has a semi-height of 1. The inlet 
Was taken to be x = 0 and the exit was at x = 5. M = .44; 
the inlet Mach number on the nozzle center-line was chosen to 
mitetd Sonte conditions in. the throat. 

igo yselGerous, were possible for Ehis inlet) condition, -mthe 
symmetric solution and the asymmetric solution; but neither 
solution was achieved by the finite element method. Although 
the solution converged for the subsonic case in three itera- 
tions, center-line Mach numbers deviated significantly from 
both one-dimensional theory and from Oswatitsch's approxima- 
tion [Fig. 5]. When the local Mach number M exceeded the 
inlet Mach number by approximately .2 (M > .64) the solution 
Was invalid. Differences at the center part of the nozzle 
are due to an essentially incorrect free stream Mach number. 
Parening the solution atx = 1.5 would improve the solution. 

A second test case was run for the supersonic section of 
the nozzle with the inlet boundary on the sonic line. The 
exit boundary was left free to float. .Here the solution was 


unstable and no meaningful results were obtained. 


42 


h - i 
} doaW Igaol to 
| | rion ony 
ramib-aws 3s 2928) Ses7 nay aot tnd s olssos oft 
¢ 
beniteb yiabawod ade Artw ¢ sisson Wonstsawed tam 
| — —~ 
a ae 
Tr ed? .1 to sigiont-tese 6 gpd ¢1l = 2-28 se0tit3 ody ov i 
oh v2 = % 38 26 3ipe S02 fae oD = xed co tas & 2s 
oy neeod> eau Sabi weenes eizsem sid) ao Seaman Haah telat od: 
| segin? sir a2, 2 naizibaes stave” obs. 
oi3 -,fofsshaod Fink eheg ta2 visti” 3 otew enok2uloe 0 
saris ud ;rotsulo 224 3¢omn yee ond bne Apiautoe ; 2 hrs 
dgpod7tA better taste ts Oriel? Shy ye sevmieoe eo mol ul 
-#t42 j iit ni geso obnoedse ong 10% bears 2 we rot be 
aos? visneoktingte tesetvgh cae AvaM anés-c0 ve 2 i) ). 
~pmit xOtags 2 worismwad nmoy2 ba Viosd teacienensh-en no dye 
oft bsbswoxns 4 tedmarti oa Lest eta asdf rie 0, 
noituide wd (ea, 2 aewy Ree OR ay 6 xd naa me a a 29: 
f 7 isnee ot Zo tale 7939083 ad? are a ji leval 
eo: 5c mestda2 eba2 a 
ittu! si? gvotaet BE fe) Com - 
20 aor: bmoeteqie aieas 
stl enti oinoe sige 


.baarteside axe 


“siTnsoy 9281S Apeois YIM Ss i[Nsoy UoTeIB9RIUL OUT], Fo uosTIedwo) - - 9IN3TY 


PLOYD % 
a GO 38:0 = FZ 0 90 GO b’O SO A 8) ine O 2-0 
C80 
Ox 
Cll 
7 
UOT}EIBZOIUT OWT] + er | 


uotyntos Apesys —— 


IOqUNN YIeW TeD07 
43 


“O0Sb°O LOQUNN YIeW FOTUT 


O'S Cb O'b ce¢ O¢ 
bt 4 pe 


uoTjiewLxoidde JUuoWeTY oTUTY + 
uotzewrxoidde yosit1eMsSC 
TeuOTSUSUT p-39uQ — 


+ 


Ft + eties + 


‘SUT T-19}U9) OTZZON - S 94Nn3BTY 


IoqunyN yoepy 


44 


: 4 ~ | . a . 


_ lett eheai haat _— 
Psi! RO ep ye ett awed ; : 


. = 
hol 

ee ea: 

Kohisntaowyqe Jnemeli oxtail + 


IX. PROGRAM MODIFICATION 


the tinicte element computer program for non-lifting air- 
mattis, as developed in Ref. [5], is separated into two parts. 
These have been designated STRANL-I and STRANL-II by Lockheed 
Corporation. STRANL-I generates a finite element mesh to be 
used as inputs to STRANL-II, which assembles the finite ele- 
ment equations, applies the boundary conditions, and solves 
the non-linear system of equations. Detailed descriptions 
and instructions for the use a6 the two programs can be found 
in Ref. [5]. Only the modifications to the above programs 


will be discussed in this section. 


A. UNSTEADY EQUATIONS 
Modifications to STRANL-II to form and solve the unsteady 
finite element equations were three-fold: 


fy the new elemental matrices SC; were calculated 


and SM; 
and assembled. 


j j 


Zpervubothe Maerices were stored on an external magnetic 
disk to be accessed and reassembled later because of 
the amount of space required to store three large 
matrices, in core memory. 
3) Maevetrective stitiness matrix SE,4 and the effective 
load vector RE;; were assembled, and the system of 
equations solved. 
Several existing subroutines in the original STRANL-II 
program were modified to assemble the damping and mass matrices. 


These include subroutines NEWK, EMTC, DERV, and EMQT. Two new 


subroutines were added to perform the other tasks. 
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B. MODIFICATIONS TO CALCULATE THE MASS AND DAMPING MATRICES 
EMTC in the STRANL-II program calculated the elemental 


stiffness matrix by numerically integrating the equation, 


Si = Q,P; dA 

Equations were added to EMTC to perform the additional numer- 
ical integrations for the mass and damping matrices. All 
three matrices were calculated at the same time. EMQT 
assembled the elemental stiffness matrices for a quadrilateral 
and trapezoidal element from the contributions of the tri- 
angular elements. Mass and damping matrices were calculated 
in the same fashion. 

Subroutine NEWK, an original subroutine in STRANL-II, 
which assembled the finite, element equations for steady flow, 
was modified to assemble Dias and Slag s A new calling arcu- 
ment, NMAT was passed to NEWK, which assembled contributions 
from the triangular, quadrilateral and trapezoidal element 


Matrices into the global. matrices S; ce and Ss 


j? 
depending on NMAT being 1, 2 or 3. 
i subroutine STORE 

Given a non-symmetric matrix stored in a banded node, 
plus the Tight hand side wector, subroutine STORE separates 
this system into two matrices and stores them on a magnetic 
disk. Figures 6 and 7 show the decomposition of a banded 
matrix into banded storage, and the further decomposition 


of this banded stored matrix to two smaller matrices by sub- 


Eeutine STORE. In these fisures, D,; L, and U represent the 
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diagonal matrix, the lower triangular matrix, and the upper 
triangular matrix respectively. HBW is the half bandwidth 
and R is the right hand side vector. 

STORE requires an additional work area one-half the 
size of the originally dimensioned matrix which is to be 
stored. 

2. | subroutine TIME 
Subroutine TIME integrates the system 
Sane: + Says + eile ge 2 = ee 
by Houbolt integration. 

TIME reassembles the three matrices which were stored 
on the magnetic disk to form the effective stiffness matrix 
and the effective load vector. Once this system of equations 
is assembled, a banded equation solver is called to yield the 


solution for 9. Prgure 8 15 a ischematic flow chart of 


Tera | 
Time; |in Fig. & when 1.=1, the lower triangular matrix and 

the diagonal of the effective stiffness matrix are formed by 
adding the appropriate contributions from the stiffness, mass 


and damping matrices. When L=2Z, the upper triangular matrix 


tS tormmed in like fashion. 


C. CONVERGING-DIVERGING NOZZLE 


1. Application of the Boundary Conditions 
Regardless of the type of problem for which a set of 


system equations have been assembled, the equations will have 


the form 
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in which Ki; is annxn matrix and X 5 and Rs are; VeEeCLors ,of 
length.n. These equations do not take into account the known 
values of xX; on the boundaries. However, for a unique solu- 
tion of the above equation, ae least one or more nodal var- 
lables must be specified and Kay must be modified to render 
Z€ non-singular. For each equation i, either xX; or Rs must 
Be specitred but it is: physically impossible to specify both 
x5 and R.. There are a number of ways to apply the boundary 
conditions to the equations and when they are applied the 
number of equations is reduced. However, it is convenient 
to leave the number of equations unchanged to avoid major 
Peseructuring of the comptiter storage.. One such method is 
described below. 


fev is the Subscript of the prescribed nodal variable, 


th h 


the k row and the k*” column of the original hes matrix are 


Seelto Zero,. K 1s See to 2 and Ry is replaced by the known 


kk 
value of Xp. Each of the n-1l remaining terms of R; is modi- 
feed by subtracting from it the value KX: This procedure 
£s repeated for all the omtaetey wales... Of Course, when the 
matrix is stored in a banded mode, the algorithm will differ 
from that tor the nxn square matrix, but the procedure is 
similar. 

Subroutine BNDRY applies the boundary conditions for 
the modified program. Setting the option parameter IOPT(4), 
fae StTRANE-I1, equal to 1 wili call BNDRY which will read the 


boundary velocities and apply them to a banded stored matrix 


by the method described above. 
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When the value of the Xp is zero, as in the non- 
lifting airfoil problem, the algorithm becomes simpler than 
the above method because there is no need to either set the 


th 


k column to zero or subtract the value Ky from the 


right hand vector. 
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In computing the flow field for either the non-lifting 
airfoil or for the Laval nozzle, the following procedures 


were followed: 


1) The desired mesh was sketched with each node assigned 
a number. 


2) Appropriate input cards based on the sketch were pre- 
pared and supplied to STRANL-I to generate the data 
ama punched cards for input to STRANL-I1. 

3) The above punched cards were supplied to STRANL-II 
with three additional cards as input parameters for 
each case, plus additional cards for the boundary 
weloci ties, iG .the nozzle, solution rs .desired. 

4) Results of the finite element calculations are printed 
after each iteration, and the converged solution is 
punched on cards for possible later use. 

A. TIME INTEGRATION TO A STEADY STATE SOLUTION 

Test cases for the integration of the unsteady transonic 
finite element equations were conducted to calculate the 
Steagy, transonic, flow over a 6% thaick.circular arc airfoil. 
These tests were made using the same airfoil, mesh, and free- 
stream Mach numbers as Chan et al. [Ret. je publashed,.. These 
conditions were chosen to provide a source for comparison of 
ene results. 


Freestream Mach numbers used in these calculations were: 


Mn = 506 Sy Gsubeni.tiucasl) 


<4 


.861 (barely critical) 


M. = .909 (supercritical). 
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Each case was treated individually with = 0 used as 
the initial guess for each case, whereas Chan et al. [Ref. 5] 
used zero for the initial guess for M, =.806 and then used 
‘the computed results as the initial guess for each subsequent 
case. 

DELL, the value for the time step, is input by a parameter 
specified in columns 41-45 of the second card following the 
title card for each case when the unsteady option (IOPT(6) =1) 
us “selected, 

1. STRANL-I Program 

a. Input 
Input cards used to generate the finite element 
mesm are listed on the next page. Cards were arrdnged an 
aecordance with Ref: [5] < in thes following order: 
Excle card 
Option card 
Element cards 
Card for the total number of nodes 
Node coordinate cards 
Card for the number of boundary nodes 
Card for the boundary nodes at infinity 
Card for the nodes on the line of symmetry 
Card for the nodes on the airfoil 
Caras. £or the slope of the airfoil. 
input cards to STRANL-I for these tests are listed 
on the next three pages. 
be Output 
Output from STRANL-1 1s am the form of printouts 
and punched cards. Printouts from STRANL-I are listed on the 


following eight pages. 
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2. STRANL-II Program 
aa SpE 
Listed on the next three pages are the input cards 
to the STANL-II program. 
be OUeput 
The output from this program is in the form of 
printouts for each iteration and punched cards for the con- 
Werped Solution. Output trom STRANU=E) tor Mo -=". 909°"1s listed 


on the following eight pages. 
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B. CONVERGING-DIVERGING NOZZLE 

Two test cases were run for the converging-diverging 
nozzle. The first case was for symmetric flow designed to 
y2eld Sonte conditions in the throat of the nozzle. The 
second case dealt with supersonic flow in the diverging sec- 
tion by starting with the sonic line as the boundary of the 
nozzle mesh. Oswatitisch's two-dimensional nozzle [Ref. 9], 
y = bs +,/2(x - 2.5)4 Wied ssiema-tEhroge height of 1 ait 
i_—eeAsseWas Sead in both cases. Boundaries for the subsonic 
mezZte were at x = 0 and x = 5. 

I. symmetric Solution 

This problem was analyzed using the mesh shown in 

Pus. 10, which consists of 126 elements and 152 nodes. In 


ene “second, card (the option card) IOPT.(4) 


1 anducates that 
men-zero boundary velocities at the imlet and the exit will 
be read and applied by subroutine BNDRY. This option requires 
Pit the number of inlet and exit boundary nodes be specified 
i columns 56-40 of the mext card. Perturbation velocities 
at the boundary nodes follow on the subsequent four cards. 
Subroutine BNDRY reads u and v respectively for the first 
boundary node and then continues reading u and v for each 
igtet and then each exit node in the.order specified on the 
appropriate card. 
7. supersonic Case--Diverging Section 
hie mesh wsed for this case us sketched in Fig. 11 and 


iipue Cards to SIRANL-11 follow on the next page. For this 
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Gase eEhe options in effect are FOPT(4) = 1 and IOPT(5) = 1 
which cause non-zero boundary velocities to be applied to 
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Ses TFODORHAMNNMOVDONAMNMNEDDONTOMOEDHONTO AODKORe 
Po = a-4 rent a4 se | a4 |} j | 
WILL TL WLU Lu Wu 
Catrall. DOM 
ANAN NS DDONTOMNMSASON FODEDAMNTODORDASMMN NOonga~ 
NANO ee eo Se eS aoe TIAN & 


eeeve tal ce =e ol el =e Ce he eee ee 
ANAN Oat 
it 

ANON OSDRONSFOOEOASMNTODODAMNMNTODONRAMWTt ac 
leielele) an ee ee ee Or ae ae te Se Oy OI eer a o°0°0 
ree -—- aH ce loon Keo = i if 
WWI Ui LWW LWILL 
Mths =O 
TOFD QDONSIMNOEDNOONTONM DSHONFOWME DAMN TOO LAL CO 
ANMN Ree ee ee aS Oe oCoorFT- 
eoeegs =e as - oe are Cn eeoeoes 
DuUrauwy OOVU0s 


NIN DDONSI OOK DDONTOMNEDSONTOWEOSEMNTOD Homo 
ovUu°0o°o a Poe Thee CE ee eb OO U0O 
| ie | _ 4 — ado | ! | 
WWW Lu WiLL Wa 
rorny OOM 0O 
NOND ORDAONTODE DAMN STOUDODHMNFODONRAMNE™ DOCTORS 
NWA ee a ee ee ee re re NNO O 


oeeve a= eH | a a - eee 8 @ 
ODNON Panam 
it 

MNT DAMN TUDO UWS MEI PO DON SAMI OWOCNTMWIM™ OW bee I aah >) 
SCo0o0 A rTor sa aa, MOOMNMODS MMOMNNNOS eal eo] eo) 
Oe er aa = a ae ~ me = 1 
Wilt tL LLU Li 
m0 40 eyo 
FNFN DOVNSFMMOMDSOUNPOMNE DAIMNSOOEDAAMNFODO U\LA CO 
00> ANDO ea hel cei NOWDMNNMOD4H N™-OMNNKNON CO--ESs 
ee ee =~ = ae es oe heed et eoertee 
Meyer) O9000e4 


ONIN DONFOMY KOSONTVOOMAAMNT VDDODAMNnFODONSFOAA eS 
SOO OQ  AOOSATHODNMNNDAIIMN COMNNDAIMNSOMNNDNNOOOOO 
ee ae Gan si ~ a4 a4 am o— 4-4 | Pater 
WW LW UI LU LL Lu Ls 
TUTTO Mat coo 
SOF MPAA MNTODODSMNMOMODON AMON DOONTMOM DASH TSW 
=O RFQ APTOS TOHODMN SND ATONOOMNNNDATMOEONNNDAMNRAOON DO 
O-exeree ee = ~—— kame) ee oe —— swt © 8 & © 8 
erasta NAO 


FNS NE DAMM DO DONS DODONS MUS DAHON P OS DAM ao 
ooo0c0- Sa a el a Sr ee re eS I A: olale) 
a eat a a4 ae — so 1 ij 

WILWLI LL WU Lu 
BWOW al ee OR 
OM OMONTODEDAMNNTODORAMNNTODONASNNM™ OCDONT Wellavee) 
ee eee to ee ee Oo NT Ore 


uve et @ a aed onl aA = vs e@ 8 © 8 @ 
~PONON OoOO004 
ete hf 
SS OOK PODER AAMNNTODODNAMMN FODONAMNHEODONT ey a 
= Se ea ere ee ee ee a es oe eo (=) (>) 
| | a o Lo laen! a4 cad a ! | 
w WwW iw Ww 
oOo Oo Co ML 2) 


DW NODSONMDDONAMNAMS DODONFOAMSAAONSFOYEDHa 9 ~~ & 
ODVODDQO LRPONNNADNSTPRONNOYNSTREOFTNVANSTYEATONOMO 


Ceee ae a ae ee a Comal —“— ¢ © © © 8 
AaO-O-R joloolo) 12) 
oO 
° cy NNSMOE DDONTMOOE DAHON TOOM AAMNFODHEAAMYW a0 

QO OFF QUODN AN MODMNUIOS MOQOMNUVOr Mimo) el ele) 

| [= = = toon) ~ a4 - et tot 

Us Lu U4 Wy tL) 

ryt ~emo 
SM MON FUDODSMYPFPUOUDOVNSBMOEVS VOONTMUINE WN DONT WO uU\UVo®O 

OWMNOWN NWVUDMANDs NODMNNOVee M-—OMNWNO Ss FTrPOomMoorr-« 

8 ¢ on coal Con) arcs La | ee a meq ee 8 8 @ 

Oonon ele Neer...) 


81 


meu ethos dus at Ahaha 
to bee . — —— 
to QU Si epee 
ee ee o-ViGeviangee 
a — 
1 ma gy) Clete Cary Fy ree [e + 7 t to ey ewe 
o> em pet lar! 1 ia sl 
a — ¥ bateland 7 
, erat) CD Ua ee 1A Ls ea OE 
- Mame eh Oe oo Sereean Dats es 
— ea 


j j all | ' TAOMNSa wm- —erhtiIeogm —eh 
_ -— - 


a += 99 Db (hte Fete ys OO Oe ea OD " 


ri > Oe pin el OO roe eA Den Ge wun ovarBe 


- 
RY On “F o> tenant eh n cole 
ro | at mn aS LO op wings 


— Pies, tome. oo ew i he 


th oa - 
— 
5 GO AUTO 9 <1teer $e @-aeermual 
ro AA SPUN 9 OY Seen EY Oa GP Rea 


ay heey — _— to vr 


ee 


MITE |) ONG t- (Pitre g ODO SMOG 
wire = Tm ee mrt rd I 
- —a r~ ~—, - wry 


her “fue e1eg on ch i aw 


~ + Re (OD eee Te Oe a 


a th, — do = 


a 
o oe Ow 
_ 


OL? PR te AY OF Oe ie hey Genesee ry 7 
n “Fur G Jie (PED OM Oy oom I CUT f 
iin Inte — oe —— = —t 
| — 
a. a ee 
SU 9D) erry RS GG Te ee 
h 7 mesa My = epithe 


— _ von ~_— be — - 


ta Oe +9 SOG DAM otha 
E 2 te ex (AMES Mo sah) Ok Lee tah 


- oe my, ~ —o tomy 
é a 
> 


7 


é . vv 
FIM mA ODD H alrien OD Ce DVI ee 


. — -_ me wh. ne in 


+ ee i A) Le TT ee hE ee GS —_— 


me Pee ews tne. 


ADSjiADed 4m gee STAD NINOE 
a pms Tl) Li re ee Ct -_ Le Gye te ee ~ 
- - — =e : 


“ ~ tT PVN ite “ray SEBS “OO 
al —w wy hist wi eve 
yt -_ 


, ay —_ 
. ‘4 ~~ wt eo op wee ie Bei mew _ 
) “sim Ler ae Rue 
sade a, 
7 _ 
7 
? oF 
7 i 
ns > = : 
- PAL 7 - 
i * = - 
D -_ _ 


- 


te q: iJ = 


ORO ROR OF OFOROFO ROCF OROFORORORORHORHON® 
ars Tk oe et Pe 
1 | | | | | ' t | foes 
US dd Od) dd 0 0 dd a us a Ws Wk 
SM FINE ADDONDOVODDOONDO-R AINE Tf Seaeano 
OTDOTORAAOVODODODOVO BHOOROTTOEOMNDNOasS 
ADO FOMNOMONODNONDNONDMOMOTONDODRHNOANTAIN 
eeeesveeserteteenrees © & 8 © 8 © © © © © © © 8 & OS Bea} mat 
AT SPT IT SP tt Sar ttre re reneord 


DOCOCOOCOC COCO COON @MWO9O 090000000 0000000 
lolelelelelololelejelelclolcloloe) clelelelelolelolelolelelelelelo)'s\._.) 
a= 
UU 0d dd 0d 0 Lo) WL 
MOIANODAYNNNNANNO OC ODDDWOQOWOADNNE MY Mt too 
DVO FTOONVUVNFTMNMOO OOS EH MMAMN Nat OO OW 
ANN D DOO NAM MFTUOIMWIY WOME OD Sa4ryrndaonmMooro 
oseer eevee ee eee 8 2 © 8 8 8 6 © 8 © © 8 © 8 8 OO Heri med 
SS SSN ANNAN NAN NN ON ON NNO NN A) 09909 09) OO PTE 


Set et a tt tt tS a Ht St St i et Set a tt et Set HO RHOOr”™ 
OOO OOOO OO0O COVA VO DO OCCOCOCCCOCCOCCOTm 
eta okt ee ce Hooch SE I ere RS IS ke A est 
UU dd dd dW a a a Lk 
DOF ODAIDYD D&S FP) OD SOLAS) DENN OD OT ODMOOO 
FAD VASO M OH ODM WY DI ODE MAME NAD aSHPOODaMroD 
RONAARDAM DDO DINANDN DNDN DADODH! DHRONMBOANN $O 
ome 8 emeweoeoeeweeeWwenseese ee @esseqeeeoeese @ e © C4 
DOAVDND NON ONDNDNDNONONDNDMODHN Hera 


VOODOO DVVV90 OO OV OVVO VUVGVVOOU 0000 VOVO0OU0UrF 
eleleleleolelelelololelelolejlele) Sele) clelelolelclolelelelelele °°?) 
ee 


UU dd Os 0 WL We Og Qa WU DL Lu Wo LW Wa Wa 
MOAMDANNNNUNNNNOODDDDDDODNUNEE EMS Trroo 
DOO FFPODNNS SOMO OONOE AMON OAM SaInNNOOMNn 
FANN DDOOIONAUNANMHSFONIN DOE ADSAaSPHPODNAMOOCTE 
eeeetkeereteeereeeeesrt e?eereeeeepaee © © © © © © & Ont 
SSS SIN NNN NN NN NN NNN NANNY OY SEP ID 


SSA SS SR St SS RS SS RR HA SS Rt StS SSO 
slololelelololololelololololelolololelolololelololololslolololelownre) 
Peek ve re ER PRE Eee ER at 

WU Oy 00 OO OO Oo Lo oo 
DANDRFIOHDODNARDOANHDIOINAODODDADNDOFTOMe4 
DPPOPAR FON ONS NTA MOLD OS ETOR SAREE 
NV PPtN PSP TRS TSP PSPS TENST TOE ONODN SO WY 
one C18 CC etree @reeset fteeeere#ef 88 Cert eo + @ @ 
Oda ay aR ae ah SR SR SR SR SR GRE GORA 


elelelelecloloelolelelelelcalolo)lelelojolelelejoe)lelelolejlelejeleleyele)) 
ielelololololololelolelelelelololololejclelelolelelololeclololololoma~ ys 


QU OU LOO OO 0 OO 2 2 oo 
NAN DO NNNNNNNNOODDDDDDDONNREEEYtOO 
DOOTPOONNETOMOCOUNOE EMM AM aeNNOOMNA 
HAAN DDOONNMMTTONNNNOOR EADS DOMMOO + 
eoeeeet®eeeesre#eeeese#e#e#eee#eee#e?e#28eee08eeeeeree @? @ @ 


St SIRI ONIN NNN NN NN NN NNN ANN 909 0909) OI PLO 


on a8 
Sono (OOo 2 oO) ooo 2S 2 2 ° Oo : aa) 
bt | | | | | ! | l ! | ! | | 

EE CES CE CEE ees) sy) 
NON WT oOo = += 2 Om NN ON YY = 
moO Ww wn _ ‘onl wy WM - OH 


N N ion N N ing) 
OOF ODODOPY OF OF OF OF OF OFoDOVORnhOO0OMON, N 
ose eoeeetre ep eevee ee teerere eee eevee eee 


WBDWONIONCNHOMONONONOMNONONOMS DO0OrOD 


DO OCOV CVV OVO VOCVOOC OVC OV OV OVOC COO OCV9OON-™= 
VOVIOD VO OVO O COSCO VOSS VVUGO COVVOO9O00000 - 


dd dd Wd Ws 0 dW De Wd LL Wt WL 
NAMDONNANNNNN NO ODDDODOWONNEH RET tryvo 


120 E28 S136 


112 


40 48 56 64 UA 80 BE 


a2 


DOO TOVUANNETMMOOOONYOFE FOO MO SANNOOFWOMmOTN 


FN OD DOO ODN MP POW OP DD St ST OOM MOO 
Sees eecovwwreeceeoe eee eeveeeecveeewvwote 
SS SS SSI NNN NNNN NNN NNN NIN NO 09 OO CONE SLL 


82 


OND+SOON 
=OpODRH 
a 


DOA MODWY 
SPONDS 
Coma 


TOVWVND-+FO 
—MPOoro 
— 


MH Dine  O 
aN POM OO 
— 


NDOTOONGD 
SNTOMm DO 
— 


aM OD ONS 
SASS ODO 
ce 


VONDTOW 
SAPO ™ DO 
| 


AIA OM OIN 
NTO DO 
ec 


DBFrOVWVNVT 
NFU DO 
— 


Mm ADMIAMEO 
NOW OO 
- 


ONOTOON 
NOM OO 


PTOWVNaATO 
NMANOOO 


ANDTOVUND 
Oo SMO ON 


uw 
tf 
NSE ODWah 
MaroOom 
ei mie @ @ 


AIMFoOoDd 


NON 


c__— = AR AR B32 


mh “Dy € -SeQm 
~*~ avproojcpe. 
aad —— } t a Gel 
Vit td LADD |e 

oD ee eee. Ge 

a 4) aa As 


~ ae 


ve tite Bd Ate Oe 4 POHOMOVOo 
_ eee teow bhee Oe eee ee eee 


ae aE ees Sse 


LOO DOCOCCOoOC oC OO 
cocoon Scoocooodooo 


yl eM Ua ay a dd A pe wee wanda 

o> alarhnterh ATS Gp CAR? GD ETI genes Wri fr 

Lb eee Cor etter eee GS “ 

’ D*~nCGss Gee arate ne Ac oom BAUM 
: ‘pe eee Meee veaeweeteanc tases 

de> em ees a iain art: 


7 oe ae am Me hn fe marin i mf 
ary OO Poole >So? cocwocouT 

wenn vaneceeree se 

AL Ae Ua A le Me? a Leg 

pon eS OG Vs ets ee ee 

“YO bre"? 4104 ee fee aie eee ae aN 

“he fe OR LATS AED er ie watosencrre 


(tee rn tp ate a oes -e* aearn 


~ pe wer Qin ek 0 ook taT * = Merete taretey ss 


ar > wi Duc OU MEU ERAIODY i 
‘oO ee 
_ 7 uray a bus 

WElPTivebe tame leeyetis tir) yee le” Pie! bb WU 

Ott > jt iY eS &@eweooe. 


aed | aw “yr = oe Pieri OO 
at SO ete ey me ors 
‘ **renee@eeeeneaee 


‘Se? SMO MmrrMmmrmy sy niataaintaingnnaintan 
ee oe A tect mee hy te 


2aqo 2 —VeCobooo Saat 


PA Re MURS, Oe? 

wd tA Hidddal tet) Apia mpeee Lee 

ine = ¢ "Os pom eo Ca laet ee 7 apap . 

> f ae § 7 it ; A * : A 
oe «@ hints ees > on ord ‘* a yi i 7 


ee 
7 mn 


' 
be . 
fm 19) A ey Ce me po eae 
v 


OoCoooVveosod COMoOeReo ©2ocoeo 
OVO.JOOO IOVS 


Lh eb bapa 7 bet 

—— * we ot wi MorwaceaT 
— wo Seheee ih om Furie Auth dl 
t - of © Bat] 4- tee AQ 4 ie G» “8 
oui ss, rea ner eae ** 


PUIG em PC CAT 


- — = — 
3 oe eC. 
; i ' ‘ ‘ ! ' 
view wm & ee i 
ewe wm 
a" a 
- 
. 


hm 


~ 
—_ 


of 
oT 


ir wy A 7 7 
Sr Se CD 
“ee tate be 


" re D6. poner rt: stiaeniaa 
— wa 5 


> —WIGOU OUVvOooVvVeVe ESSS= 2c 
 - SOU CPE ooweogo™ s 


- we oe - x . 

wy i ann a aad an ee Wh id sascaaahaidalga UL 
Oe & Wer way th Hat to & 2 ; 

a or Ra Doe at ind ne 

2 ramet # BOM Ware ye ese ey 


oe 


Ae SMM 
? 


NS QOvVv0O 
a vot ht 
WIL UW Wu 
Toto 
~o oror 
NS NOND 
=o eee 
SNRIN 

rd 
ON MANN 
NS coco 
a4 mt wt 
USUI LY LL 
oror~ 
na ~OrmoO 
NS OVO. OY 
=—=— ee @ 8 
de el ed 
Hot lat 
To NAN 
NS Co0°o 
—4q ete De a 
@ WLU Wu Lu 
oe) TON 
mo oO oOmMmomMm 
NO oror 
=——4 e ee & 
1} DNON 

= | 
NW MITA 
NOH Ww leleleie) 
aa ee ie ie 
| LULL) Lu LL 
N MLA 
-—~- OO SANNA 
NY 2 TAT Y 
-— eevee? 
& MAM OH 
= feat 


83 


« - ~- "7 
é ° 


: J 


mire —<«,* lie 
WwW) la" 
Shi pity * 


INPUTS TO-STRANEoTT 


UNIFORM McSH 


»SONIC LINE INPUT, 


STEACY SUPERSCNIC FLOW 


MAN TOEPRONMMIEDOSEOH TO BHOKHOROHOROROHOHORHORO 
rare MNONINOASEE OM MOHON CDOVDOVDOODVDOCDOOOCD0O000CO0O 
| t | k t f | ' | t 
LIL UU Ud Wd db Le 
ano OM FODAIVNGM DOOFMOM MDS O 
NO DAODPANLPUPDONMWN SAHNONANNMYDXOQAMONADOUWO 
NM DNONMISOASMFH MONONA NOMOFOWDDOORMNNE MO STN 
ee eoeeesrses?eg 8 eeene @#e 8 © @ @ @ 
on N DY ee Dw oe ew a oo nem lat let liete Net Nel. 


WD DOAMNPODDRANTFO™M NGA NAeAesORnOOOOCOCOD0COCOO 

NOSOaT NONDONMNAST DODDDOODOOCOOCOOOCO000CO0OO 
DG UG) Dot if 

WO WU UW Li) Wa bo Wu 

THYNOMOTDOMNUNTFOVUNDOVO 

OD DANKMYUTASHAONMNOD BAPONANIM™OTOAODDDNANDOO 

OO NOSIS MNOMONNAT AMNDDA™ OM AOTNONDONRINWY 

ee - eoeoeteesteeseeteeeeeeeteeetee @ 

oto) LANM DOAN DOA DSS HSS NUNNN 


! 

MAM DOAN TOE DONO ODO SS HHS ASSO Hn OnFOCHROFO 
ord NW SNe FM MOMVNNATOOOCOODCOCCOOCOOCOOOO0OoO 
| aes att ith sik ! | | ' 
WIL dd Lo 
ru WO DONO MOMS OM ONQUNATM-PDAIAM 
IFN OSDVOSAMNNOOQODKRANLTNOS DE NHONHNDD DE HNHOOCOUNMNHOD 
OO NWSI TE ONONONDNASHE DODOBDEOONS TODO DHONMAN 
ee eeeteethrteeeoeostee*e#eet oe @ © @ @ 
uwVOV NONDNOMDIMNDOSOMI MS fetrSs 


MmRKON MEE DORM TVUPMP NAN SS SS S434 O0O OODVD0O0VO900 
aTraeT NONONNAT MKCOCOCOOCOCOOCOCOCO0O000°0 
be te Ba 
UU Od dL Lt 
DOOM ANMNVDANMODOAPNNODOO 
DOHMHONFODDANSFONE DONDAMNAMNMNHOIMNNDSHOAMNOCOO MO 
SO BMI F HNONONOSF SEOTNNNTAMNAOTMHOODDANNNN Oo 
ee oeertee tt epeerteete eet ee ¢ @ 
oo DAM SOS OP BSS SS sR AINANS 


| 

DMT DOSIMNCDDHNTOAE DON SAAS As SASSO FOOMNN 
rt esta jee inte 2S) at ioe Utes ne Jc TAL Vis [ee He \ 
wilu CU Lb Li 
Tr DNO~STFHPODNOEONEYETORRME 
DD PHONONS DOFIMSFODDAMNDNARDON ANS COCR AONE SErst 
Sa OSI TATE MON ONO SEE ME TAIT NTMOMNOODMTFODDTFOORO WoO 
°° oeoeeeteenrte#et eee e#eeee 8 6 


aie AM Ay a aN HO RHOHONANG 


| 
4 DOANTOMDONMNNDODORFMN AN SARAHB OO DOO OOCDOCOCOCO URDU 
tT NONANMNSAS M~NOOCDOCWOOOCOOOQOCOO0O0C0O0O00000 TN 
Vole late t 
QU OU Od 0 0 0 od Od 0 oo Lo 
MNSOOPFIANRORAMODOSNOO 
DONMNMSDOSMOHTOP DANSON LADODONSES HOOTONMNAHOONHaO 
CSO SBT NONONANS SE NENHOOMENMOOD BOSOM ODANNWWM Sot 
ee eee? eee2#:e2e6egoeeeeteereetet ee? &® 8 @ @ 
ooo SOMO OM Ae sets NUNN 
° | 
Sa NOONTAEDOSMPODRANS BP BA BA BA A ae BAB KR St AEMO 
OD SAT eMONONMNATE MEMO CO 9 4 ol) (iis) a Se & 2 m+ 
! = | | | | } ! 


E~ 
E 
EF 


os CE ee es Ww wi 
0 OA N Ps, 7 f= a “a Oo 
BAGS OS NONmAR Sem Se (8S ONY ton) oOo O Oo AMNON 
oOnwo Som odindines ari mOCNoL ooo sos oro nO GOdGm NO 
*¢ eeoeteeeresnetersege??  ® &@ @ 8 @ 
“or ODNONONODCODOOOF-OrF0ODOO 
Oo 


° OSMNDO DD ASNTOAM DONMMOAMN Sessa aOODCOCOUOCOOCOOCOCOCOONME Tr 
DW MONWOSAST MVOWNNYOO OV DD OVO VDGO9D9D0C C0000 main 
Det Stat 
WU 0 dd Wa Lb by WL 
WNOFOTINAINNODNEONWINOO 


La MNTOMDONMWDODO AM FPOMN SO DY) VU0O Tr sTDO OR TMONOU MNO 
OoOovVW MOUNWIN Wit COR OOD) tS DP TW SM SPM ODOANNYYY = 
ee eeeoveer vee eevee tr eeve ee 
oo SSN NON DO Sass SS ANNNAN 


84 


HTS OS OSOWOMOMG : 


te tee thoes ee OS ES PIC Py belle 


i) ee ae ee 


O° Doc OPN 
j ¢ —_ 
aA Ad Ad Gl i ee wpels r Aas 
- ~ +4, 7 fat?“ tr eee 
£ UL ar dil ap  — pe es 
5 (2 Soors som stamina tien 
* - * *fFr te @Ceeveet ee - a> ne 
oe Te Sm 7 oop 
ot nm bee et hgh A ee a J 
BODO Soeodeeco i RNA DING & i — joey 
’ | 4 ¥9 af LI 2 at s 
ijt ad ASU Ly) _ S : 
“ts "OD YPOnmay ia] - 
oO SAS =O feito Dee he baited _* 7 
1 TRE is Ee ATO Behe mye: CD” 
*% > ee ® ee | @ s/ . , 
my Gerri a 
> i - > teh ODS: HRSG AE Pare VS _ 
od soo CoS PO AIR Te OY he Pye tn le Peeere ~~ 
i ifreTre® e' | 
is Ce 
: ‘se “Ru |42ore> & i ToA4y 
veo NaPA ea haan Ler en OG UEDA i = oy 
‘ t Se oS 06 24 = tO Ore Fae om 2 
. . * « s Pet? eetege«é os ie 
- qa ene pen we fhe ice 4 
a = 
t Phe _ pean ieee pth ceteeh Nie eh odbt deadi) hep aks ae) s. 
; 2aKSOoAh SKGTHOWOM peta 2 
ee oe eee 0 oa 
bas a ids aid la Cl 
i EDT og COV a = 
*] te an — UD een GOD A ti CLO bike mtge © VS . $s 
ror Vite Sere 9S to MANO Ge tee Alaa oo -£ 
» © ' » > 6 FF 68.2 & es 
De ton FT ee ‘oe 
= 
layne es ee conte es ey pte a Nee A CG Cee 
~ i 2~IQCogo- OAM OTe wales o J 
‘4 repigra _— t *] 
? NSPIPYEROIP Tees eel val tee a as 
1 OPA SB eet > — 
Dr? LS pte 1D > ee Cod ) =>" TO MOAT PIO Bag oo 3 
DCG t- "Dc OTP a — i 
’ * : «a ® .> « &'e*#tte*® 
- p= Oto — pes hme ote } ott r 
- 
one: nm or Oy LP De FPS OD et OMI? CIT om ~ 
OHOoSOoasS SACO a mam preite At oy apes a, 
>t 4.8) 98 J 
(abst es (2) artes CPL is 
Ab 4 > @ owes 
c (9 th ie + So ae 
—_ 


wn > Ow 
& >1 ne Berens wow 7 
mi hth + Te sane PS 


4 sila’ 


re et Be ee ee 8 @& ss € 7. a oe 
pa be tamer Soy ol PED ere Oe . —_ 
1 ; — 
— eb Sk Ml : 
O Moat 3 


> =-— mn t+ Mim Oma 
“VE 62 wwe AE a Cppabniedee 
. es @ ¢ & @& & #8) 8.8 


~S& Dires pe hens ert COD 

146 PRA gaits 
1iwe 1} 

Fee ayr A sere te ++ 


e * ” 
Steete tees a So Bis emmeiia 


UY, Wy Bone TSE EM a temo ane EY ‘> 8 


*ereaewvrer PRE @ 4 


Ea 


We SLTE-O02 1.289E-01 1,.786E=-01 2.215E-01 2.576E-O] 2.597SE"Ot Sehlie-0l 


MOS Muy 
Mey ITSO 
se 8 
mit) 


85 


ae 


WHE 


86 


Figure 10 - Finite Element Mesh for the Converging-Diverging Nozzle. 
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